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SUMMARY 

M. 

A  six  degree  of  freedom  flight  dynamic  model  incorporating  an  aerodynamic  database 
derived  from  DSTO  wind  tunnel  data  developed  for  the  Mk-HI  store.  The  model  can  he 

linked  with  aircraft  flight  dynamic  models,  and  forms  part  of  a  general  progranmie  to  model 
aircraft  carriage  and  delivery,  weapon  release  and  store  ballistics.  A  rudimentary  ejection 
model  was  also  incorporated  to  represent  the  forces  and  moments  experienced  by  the  store  at 
the  point  of  release.  The  tnodel  is  valid  for  free  stream  flight  only,  and  ignores  the  effects  of 
aircraft-store  interaction. 
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Store  reference  diiimetor 
Anglo  of  sidcslij) 
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1  INTRODUCTION 


Uiuler  Air  l  orn-  |{.■s.■ar(  Il  lUquiroiiKUit  AFRR  7/!}0  •F-lllCVRFlll-C  Modelling’,  DSTO  is 
providing  atipport  to  llio  ItAAF  in  the  verification  of  store  release  and  delivery  accuracies  of 
the  weapons  sysleiiis  in  the  F-1I1(J  Avionics  Update  Programme. 

The  acciirjuy  of  sluie  delivery  depends  on  a  number  of  factors,  including  aircraft  release 
conditions,  .store  ejector  characteri.stic.s,  aerodynamic  release  disturbances,  store  aerodynamic 
characteri.stics  and  atmosplieric  disturbances.  The  calculation  of  release  point  should  ideally 
account  for  the.se  error.s. 

As  part  of  the  proce.ss  of  store  delivery  prediction  and  verification  DSTO  is  developing  dynamic 
models  which  can  simulate  the  store  carriage  and  release,  and  store  ballistic  behavior  of 
unguided  weapons.  Information  for  use  in  these  models  is  being  derived  from  available  data 
rci)orts,  and  from  wind  tunnel  and  flight  test  measurements. 

In  this  report  thi‘  <levelo|iment  of  a  flight  dynamic  model  for  the  Mk-82  store  which  can  be 
linked  with  models  of  aircraft  flight  behavior  is  described.  The  model  is  valid  for  free-stream 
flight  only,  i>ut  is  structured  to  enable  the  effects  of  aerodynamic  release  disturbances,  detailed 
ejector  charact<'ristics,  and  atmosplieric  disturbances  to  be  incorporated. 

2  FLIGHT  DYNAMIC  MODEL 


The  flight  dynamic  model  of  the  Mk-82  store  was  developed  using  a  standard  Six  Degrees  of 
fVeedom  model  of  rigid  liody  motion  developed  in  Air  Path  Axes  (Reference  1).  The  model 
was  written  using  the  Advanced  Continuous  Simulation  Language  (ACSL,  Reference  2)  at 
ARL  Melbourne.  The  program  uses  a  fourth  order  Runge-Kutta  integration  method,  with 
quaternican  jiarameters  for  the  calculation  of  aircraft  attitude,  and  direction  cosines  for  the 
gravity  vectors.  The  u.se  of  AC.SL  enables  a  range  of  dynamic  system  analysis  tools  to  be 
used  during  model  dcvelo])iuent  and  application.  The  information  used  in  the  simulation  of 
the  Mk-82  store  is  listed  below. 

1.  The  physical  dimensions,  mass  and  inertia  properties  of  the  Mk-82  store  were  necessary 
for  body  axes  force  and  moment  calculations  (see  Appendix  A). 

2.  Body  axes  forces  and  moments  were  calculated  in  subroutine  AER02,  which  retrieves 
the  appropriate  aerodynamic  cocfRcients  from  the  wind  tunnel  database  using  database 
access  routines  which  convert  the  aerodynamic  forces  from  missile  to  airpath  axes  (see 
Appendix  D). 

3.  A  database  of  aerodynamic  coefHcicnts  was  incorporated  in  the  form  of  data  tables  with 
independent  variables  Mach  No,  angle  of  attack,  and  roll  angle.  The  original  database, 
obtained  from  wind  tunnel  tests,  covered  only  positive  angles  of  attack  and  contained 
non-zero  aerodynamic  values  at  zero  angle  of  attack  and  roll  angle  (see  Appendix  C). 

4.  Unlike  an  aircraft,  the  Mk-82  store  cannot  be  trimmed  using  movable  control  surface 
deflections.  The  equilibrium  condition  of  the  modelled  store  is  achieved  by  the  addition 
of  incremental  body  axes  forces  and  moments,  and  the  use  of  a  non  linear  equation  min¬ 
imisation  technique.  This  procedure  is  initiated  by  keying  the  command  ‘eigen’.  With 
the  store  in  equilibrium,  standard  Unear  dynamic  analysis  techniques  available  within 
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ACSL  ran  lie  iisotl  to  uiuily.sn  llio  linear  dynainir.  behavior,  and  provide  information 
aboiil  tlio  imnlcs  DfoKrillalion  and  other  tnotions.  See  Appendix  D  Example  Simulation 
Run,  .step  12,  for  an  explanation  of  the  eotninand  ‘eigen’. 

5.  Store  ejection  was  siiiiulatcd  by  the  addition  of  the  two  subroutines,  EJECT  and 
EJI’(n'2,  which  provide  a  rudimentary  representation  of  store  ejection  forces  from  the 
outboard  and  itiboard  .stations  of  a  niHJ-33/A  Vertical  Ejector  Rack  (see  Reference  3, 
also  .see  Appendix  H).  Ejection  is  achieved  by  applying  fixed  ejection  forces  based  on 
measured  data  to  the  store  motlel  for  0.12  seconds. 

G.  The  subroutine  MAKEFILE  was  included  to  generate  an  ASCII  file  of  the  store  trajec¬ 
tory  position,  velocity,  and  orientation  for  use  in  a  graphic  display  of  the  behavior  of  the 
store  flight  behavior,  using  a  program  developed  for  a  Silicon  Graphics  Workstation. 

A  more  detailed  <Iiscussion  of  the  modifications  and  additions  made  to  the  aerodynamic 
databas.?  is  pre.sented  in  Apitendix  R. 


3  THE  MK-82  AERODYNAMIC  DATABASE 

3.1  Database  Source 

The  Aerodynamic  Database  was  obtained  from  two  sources  of  experimental  data,  viz  the 
DSTO  wind  tunnels  at  Salisbury,  South  Australia  and  Melbourne,  Victoria.  The  Melbourne 
data  were  mcasurcti  in  the  transonic  wind  tunnel,  and  covered  the  Mach  number  range  of 
(approximately)  0.9  to  1 .2  ,  with  angle  of  attack  (o)  0  to  30**,  and  roll  angle  {<j>)  -45“  to 
-H45“.  The  tests  wore  performed  using  a  standard  pitch/roll  test  rig  with  sting  mounted  strain 
gauges  to  me<asuro  the  forces  and  moments.  Note,  with  sting  mounted  models,  combinations 
of  angle  of  attcack  and  roll  angle  arc  used  to  provide  the  required  ranges  of  angle  of  attack  and 
sideslip  angle.  Dccausc  the  sting  is  attached  to  the  movable  test  model,  the  resulting  force 
and  moment  values  arc  measured  in  body  axes.  However,  the  experimental  results  of  both 
Melbourne  and  Salisbury  tunnels  were  resolved  into  missile  axes  and  incorporated  into  a  new 
database  with  a  difTcrent  format. 

Body  axes  and  missile  axes  are  defined  as  foUows: 

•  Body  axes  are  a  set  of  orthogonal  axes  that  are  fixed  to  the  body,  with  the  x,  y,  z  axes 
for  axi-symmctric  bodies  aligned  with  the  principal  moments  of  inertia  of  the  body. 

•  Missile  axes  are  similar  to  body  axes,  but  with  only  the  x-axis  aligned  with  the  corre¬ 
sponding  body  axes  moment  of  inertia,  l„.  The  y-axis  remains  in  the  inertial  x-y  plane 
(horizontal)  regardless  of  the  motion  of  the  body,  such  that  body  axes  and  missile  axes 
differ  by  a  rotation  of  ^  about  their  common  x-axis  only.  Figure  1  shows  the  difference 
between  body  and  missile  axes. 


3.2  Original  Format  of  Database 

The  database  file  contains  the  following  data  in  17  ASCII  columns. 
1.  Data  point  number 
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2.  Iiu'iilciiro  .s(‘t  iiuiithcr 
Roll  atiffh'  s<'l  iiiiimIht 
I.  Mat'll  si't  niinilii'r 

5.  Aiiglt!  of  iiiri<lt'iiri>,  ()“'  to  +,‘J()‘'  (2"  incrcniciits) 

(i.  Storo  roll  aiiKl<‘,  — iri"  to  +45"  (7.5"  incrcmciitti) 

7.  M.icli  niiinbcr,  O.-l  to  1.2 

H.  Cx  -  noil  tiiiiiensioiial  fort'o  toefRciciit,  x  direction 

!).  Cy  -  noil  (litiR-iitiioiial  force  coefficient,  y  direction 

10.  Cz  -  non  diiiii'insioiial  force  coefficient,  z  direction 

11.  C'l  -  rolling  nioiiient  coefficient 

12.  C,u  -  pitching  nioinent  coefficient 
i;i.  C'„  -  yawing  nniiiient  coeflicieiit 
Id.  C'(,,  -  rolling  inuinent  derivative 

-  pitcliing  and  yawing  inoiiicnt  derivatives  (a  single  value  is  provided  for  both 
variables:  see  section  .'1.3) 

10-  -  yawing  nniinent  derivative  due  to  rate  of  roll* 

17.  CVp  -  force  derivative  due  to  rate  of  roll* 

•  Since  the  contribution  to  resultant  force  and  moment  provided  by  the  terms  Cn^ 
and  CVp  is  not  significant  and  the  quality  of  the  data  for  these  terms  is  not  known, 
these  variables  are  not  used  in  this  simulation. 


3.3  Assumption 

Close  examination  of  the  database  shows  that  values  of  the  variable  Cm,  lor  a  =  a,0  =  b  are 
identical  to  those  of  C,,,  hr  a  =  b,/3  =  a.  This  is  not  strictly  correct,  but  if  the  approximation 
that  a,  =  0(  is  used,  the  values  in  the  case  above  are  valid.  The  following  simple  comparison 
test  demonstrates  the  insignificance  of  the  error  in  the  assumption  over  the  range  covered  by 
the  database: 
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u,v,w  are  tlio  roiniHinonls  of  llio  air  i)ath  velocity  vector  V,  along  the  body  axes  x,  y,  z. 
Figure  2  sliow.s  ilic  difrereiice  between  0(  and  a,,  and  0,  and  0t. 

Notice  tliat  the  worst  case  error  is  for  the  case  where  <»  =  30“,/?  =  1“,  OT  0  =  30“,ft  =  1“, 
where  the  associated  error  is  approximately  15X  .  For  the  remaining  tabulated  values,  the 
error  is  less  than  0.5'/  ,  indicating  that  the  approximation  rtt  =  a,  and  0,  =  0t  is  reasonable 
over  the  range  covered  by  the  database. 


3.4  Revised  Database  Format 

A  set  of  standard  dataljase  handling  routines  dcvelojred  by  Air  Operations  Division  at  ARL 
for  use  in  (light  dynamic  simulations  is  used  for  accessing  the  aerodynamic  data. 

The  original  datab;isc  was  changed  into  the  format  required  by  the  standard  database  extrac* 
tion  routines.  Tlie  now  database  has  the  following  format: 


title  mk82cxdatabasc 


thruput  cx 
thruput  cy 
thruput  cz 
thruput  cl 
thruput  cm 
thruput  cn 


(machno.phi 

(machno.phl 

(machno.phi 

(machno.phi 

(machno.phi 

(machno.phi 


thruput  dp  (machno.phi 
thruput  cmq  (machno.phi 


.theta) 

.theta) 

.theta) 

.theta) 

.theta) 

.theta) 

.theta) 

.theta) 


machno  varbpt  0.4.0.5.0.6.0.7.0.8.0.85.0.9,0.95,1.0.1.05,1.1,1.15.1.2/ 
phi  conbpt  -45,0,45.0,7.5/ 
theta  conbpt  0.0,30.0,2.0/ 

* 

CX  POINTS 

-0.1556,  -0.1552,  -0.1524,  -0.1476,  -0.1472, 

-0.1440,  -0.1394,  -0.1327,  -0.1264,  -0.1170,  etc. 


The  database  subroutines  will  calculate  the  longitudinal  and  lateral  force  and  moment  coef¬ 
ficients  corresponding  to  the  flight  conditions  of  Mach  No,  roll  angle  <f>,  and  pitch  attitude 
9. 
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I  hfi  tlirupiit  .statc'inoiit!)  notify  the  siibroutiiieii  that  each  force  or  moment  coefficient  is 
(letermiiu'il  Irom  the  iinh'peiHlont  variables  Mach  No  (abbreviated  tonachno  in  the  program), 
phi,  cam!  theta  inlormation,  where  M,  S  arc  Ccalculatcd  in  the  main  ACSL  program  or 
subroutines  ami  are  passed  to  the  database  access  routines.  The  force  and  moment  coefficients 
arc  called  functions. 

The  vnrbt  statemeiil  notifies  the  subroutine  of  all  the  values  of  machno,  the  values  not  being 
spaced  ccpially  apart.  The  conbpt  statement  allows  the  lower  and  upper  bounds  of  the 
variables  phi  and  theta  to  be  stored  with  a  constant  interval,  thus  using  memory  economically. 

The  values  of  the  database  funcliotis  at  the  discrete  values  of  the  variables  machno,  phi, 
and  theta  are  called  breakpoints.  The  data  arc  listed  until  all  breakpoint  combinations 
have  been  used.  Tiie  ordering  of  the  breakpoints  is  such  that  the  function  depending  on  the 
last  variable  varies  first.  For  example,  if  a  function  Cc  has  two  variables  which  have  two 
breakpoints  each,  llie  correct  order  for  the  data  to  appear  is  ; 

thus  cycling  through  all  of  the  possible  combinations. 

The  original  database  was  provided  in  an  ASCII  format  file.  The  database  modification  was 
performed  on  an  HIM  l’S/'2  model  70  personal  computer.  A  listing  of  the  program  used  for 
the  conversion  is  presented  in  Appendix  C. 


3.5  Database  Dias 

The  aerodynamic  data  for  the  Mk-82  store  as  noted  in  Section  3.2  covers  only  positive  angles 
of  attack  up  to  .'lO"  ami  roll  angles  of  -45®  to  +45®.  For  the  flight  dynamic  model,  it  is 
necessary  to  assume  that  tlie  data  are  symmetrical  about  0®  angle  of  attack  so  that  positive 
and  negative  changes  in  botli  angle  of  attack  and  sideslip  can  be  simulated. 

For  a  symmetrical  aerodyii.amic  sliai>e  at  zero  angle  of  incidence  and  zero  angle  of  roll,  it 
would  be  exjiceted  that  (litching  and  yawing  moments,  together  with  vertical  and  side  forces, 
should  be  zero. 

Inspection  of  the  databa.se  showed  that  non-zero  values  of  pitching  and  yawing  moments, 
as  well  as  vertical  and  side  forces,  were  measured  at  zero  incidence  and  roll.  The  non-zero 
terms  became  evident  when  the  simulation  was  run,  resulting  in  non-linear  behaviour  at  low 
angles  of  incidence.  Figure  3  shows  the  limit  cycle  oscillations  of  a  due  to  the  non-zero  bias 
in  the  database  with  tlie  model  performing  pure  pitching  motion.  Figure  4  shows  the  limit 
cycle  oscillations  of  o  when  all  the  degrees  of  freedom  are  included  in  the  model.  To  ensure 
continuity  of  the  forces  through  zero  angle  of  attack  it  is  necessary  to  remove  these  biases. 

The  unmodified  database  caused  incorrect  eigenvalues  to  be  calculated  by  the  ACSL  analysis 
capability  and  tlicrcfore  c<au8ed  unpredictable  stability  characteristics  at  low  angles  of  inci¬ 
dence.  Removal  of  the  bias  by  linear  subtraction  caused  the  amplitude  of  the  limit  cycle  to 
approach  zero,  and  resulted  in  acceptable  behavior,  as  illustrated  in  Figures  5  and  6,  thus 
providing  an  acceptable  solution.  Figure  7  shows  a  plot  of  versus  tf,  clearly  showing  the 
database  bias. 

It  is  recommended  that  tlie  current  database  be  replaced  by  a  more  accurate  database  when 
one  becomes  available.  Furtlier  work  is  under  way  to  improve  the  quality  of  the  aerodynamic 
database. 


4  PROGRAM  STRUCTURE 


I  lio  Mk-.S'J  .siiiiiil.tlioii  |no<fi;im  njiisisls  uf  a  main  ACSL  program  wliich  contains  the  equations 
of  motion  ami  a  sot  of  FOirrifAN  suhroutiiics  which  determine  the  aerodynamic  forces  and 
moments.  Various  additional  subroutines  provide  necessary  initial  conditions  and  calculate 
informatiiiii  for  graphical  ontjnit.  Figure  8  shows  a  (low  chart  of  the  program,  while  Figure  9 
sliows  tlie  block  diagram  of  the  air  path  axes  system  for  the  store  motion. 


5  VERIFICATION  OF  FLIGHT  DYNAMIC  MODEL 


In  order  to  determine  the  v.didity  of  the  six  degree  of  freedom  model,  a  series  of  tests  was 
performed  on  <  tuisti  ained  version.s  of  the  model.  The  tests  involved  separating  the  lateral  and 
longitudinal  modes  and  analysing  them  individually.  All  tests  utilised  the  revised  database. 


5.1  The  Lungitiiditial  Model 

The  first  lest  involved  constraining  the  six  degree  of  freedom  model  to  a  three  degree  of 
freedom  longitmlinal  model  by  removing  the  linear  side  force  components,  gravity,  and  the 
roll  and  yaw  moment  c<jmponenls  of  the  full  model.  The  reduced  model  could  now  move  only 
in  the  x-z  plane. 

After  trimming  the  store  at  <r  =  0.0",  an  eigen  analysis  w:»s  performed  to  identify  the  dynamic 
modes  of  the  K'ducod  model.  The  eigen  values  were  : 

1.  -0.00H7l.'(i!»!i 

2.  -0.2128  ±  I  .JMSOj. 


The  first  mode  is  a  stable  transient  mode  which  corresponds  to  the  drag  force  acting  along 
the  x-axis. 

The  second  mode  is  a  stable  oscillatory  mode  which  corresponds  to  an  undamped  natural 
frequency  of  u;  =  4.95  rad/s,  and  a  damping  ratio  of  C  =  0.04.  This  mode  is  characterised 
by  a  porpoise  type  of  motion,  with  pitching  and  vertical  translation  oscillations  occurring 
approximately  90°  out  of  phase. 

To  illustrate  the  oscillatory  nature  of  the  store  motion,  a  simulation  run  was  calculated  with 
the  initial  comlition  of  a  =  5°.  Figure  10  shows  a  plot  of  a,  while  Figure  11  shows  the  eigen 
vectors  relating  to  the  oscillatory  eigen  value. 

The  transient  mode  is  illustrated  as  a  drag  force  because  the  eigen  vectors  indicate 
that  velocity  decreases  as  x  increases. 

The  oscillatory  mode  shows  a  simultaneous  pitch  oscillation  and  translational  ver¬ 
tical  oscillation,  with  0,  z,  and  a  lagging  pitch  rate  q  by  approximately  90'*. 

Note  that  /f  is  always  equal  to  zero  in  this  case  because  the  store  may  not  move 
laterally. 
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The  Mk-82  store  is  iiuHlcIle<l  assuming  identical  aerodynamic  forces  in  pitch  and  yaw  as  noted 
in  Sections  3.2  and  3.5  and  so  will  have  identical  lateral  and  longitudinal  dynamic  properties. 
The  results  of  the  dynamic  analysis  are  consistent  with  these  assumptions.  Further  testing  of 
the  dynamic  behaviour  including  the  effect  of  rolling  moments,  inertia  coupling,  and  non-linear 
aerodynamic  effects  should  be  undertaken  to  verify  all  aspects  of  the  six  degree  of  freedom 
flight  dynamic  model. 

Figure  12  shows  a  plot  of  fi  for  a  simulation  run  with  initi^  condition  0  =  5",  and  Figure  13 
sliows  the  eigen  vectors  corresponding  to  the  oscillatory  mode,  along  with  a  representation 
of  the  associated  motion.  This  mode  is  characterised  by  a  ‘fish  tailing’  type  of  motion,  with 
yawing  and  horizontal  translations  occurring  approximately  00*’  out  of  phase. 

Figure  14  shows  a  jrlot  of  a  versus  time  for  a  simulation  run  of  the  full  six  degree  of  freedom 
model  released  at  an  altitude  of  1000  feet  and  initial  angle  of  attack  of  q  =  5.0®.  The  change 
in  oscillatory  behaviour  at  time  t=3.5  seconds  occurs  when  the  roll  rate  coincides  with  the 
pitch  and  yaw  frequencies  and  demonstrates  the  effects  of  roll  coupling. 

‘Roll  coupling’  is  the  term  used  to  describe  instability  caused  by  the  dependence  of  pitching 
moment  on  both  pitch  and  roll  rate,  and  the  dependence  of  yawing  moment  on  both  yaw  and 
pitch  rate. 

These  preliminary  tests  indicate  that  the  six  degree  of  freedom  model  behaves  as  expected. 
However,  it  is  recommended  that  model  verification  and  validation  tests  are  performed  as  fur¬ 
ther  data  from  wind  tunnel,  computational  fluid  dynamics  and  flight  tests  becomes  available. 


6  CONCLUSIONS 

A  six  degree  of  freedom  flight  dynamic  model  has  been  developed  for  the  Mk-82  store.  The 
model  includes  a  linear  analysis  capability  which  is  useful  in  the  understanding  of  the  store 
motion.  Preliminary  tests  indicate  that  the  model  is  valid,  but  further  testing  is  recommended. 

The  database  used  in  the  model  was  compiled  from  wind  tunnel  test  data  obtained  from  ARL 
Salisbury  and  ARL  Melbourne  wind  tunnel  tests.  Bias  terms  have  been  removed  from  the 
data  for  use  in  the  six  degree  of  freedom  simulation  application.  It  is  suggested  that  the 
current  database  be  rcidaccd  with  a  more  accurate  database  when  one  becomes  available. 
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Figure  1;  Body  and  missile  axes  systems. 


Figure  2;  The  difference  between  o,.  and  p,. 
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Figure  3:  Result  with  gravity  and  linear  forces  set  to  zero,  pure  pitching  moment. 


II 


•  • 


•  • 

Figure  5;  Result  with  gravity  and  linear  forces  set  to  zero,  modiHcd  database,  pure  pitching 
moment. 
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Appendix  A  -  The  Geometric  Properties  of  the  Mk-82  Store 

The  geometric  and  mass  data  used  in  the  simulation  of  the  Mk-82  are  as  follows; 

mass  =  kg 
Izt  =  2.35  kg 
lyy  =  56.7-1  kg  vi^ 

=  56.74  kg 

CG  from  nose  =  1.00  m 
reference  <liamctcr  =  0.273  m 

The  physical  <limensions  of  tlic  store  arc  represented  in  Figure  15  (Reference  5). 


Appendix  B  -  Additional  Subroutines  Developed  for  Mk-82 
Simulation 

Subroutine  AER02 

The  ABI102  subioutiiio  is  culled  from  the  ACSL  simulation  program,  with  inputs  of  the 
angular  rates  p,  (|,  r,  and  Mach  number.  The  subroutine  also  has  access  via  common  blocks 
to  the  angle  of  attack  and  sideslip  variables,  Alphall  and  Detail.  These  variables  arc  passed  to 
a  subordinate  routine  BULERANGLES  which  calculates  the  equivalent  Euler  attitude  angles 
of  pilch  (Theta)  and  roll  (Phi)  used  hy  the  database.  The  appropriate  force  and  moment 
coefficients  arc  retrieved  from  the  database  using  database  routines,  which  have  Machno,  Phi, 
and  Theta  as  inputs.  The  cocinciciits  arc  then  converted  from  missile  axes  to  body  axes  via 
a  rotation  of  Plii  about  the  X  -  axis.  Total  body  axes  forces  and  moments  arc  calculated  as 
follows; 


XAR  =  C,..QliAllS 

YAll  ^  C^.QnAliS 

ZAR  =  C..QnAliS 

LAR  =  (C,  -f  £^').g/},Uf.!)’.b 

MAR  =■■  (C',„  -h  ‘■^.^^).QnAltSJ, 

NAR  =  (C„  -I-  ^■).QnAltS.b  , 


where  C'(^,  and  are  the  damping  derivatives  in  roll,  pitch  and  yaw,  QBAIIS  is  the 
dynamic  pressure  nuiltiplicd  l>y  the  reference  area  (^^),  V  is  the  true  air  speed  of  the  store 
at  a  particular  instant,  and  b  is  the  reference  iliaiuctcr  of  the  store,  0.273  in. 

This  snbronlinc  lias  only  one  subordinate  subroutine,  called  BULERANGLES. 
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Subroutines  EJECT  niicl  EJECT2  • 

The  ACSL  subroutine  EJECT  provides  a  simple  representation  of  store  ejection  from  the 

outboard  station  of  the  Vertical  Ejector  Rack.  The  Subroutine  ran  be  called  on  line  from 

the  ACSL  prompt  during  the  course  of  a  simulation.  The  ejection  of  the  store  is  simulated 

by  setting  the  increment.al  forces  and  inoiiicnts  used  in  trimming  the  store  to  representative  ^ 

values  for  a  spcciiied  duration  to  provide  the  correct  impulse.  The  values  used  in  EJECT  arc; 

DZAR  =  7.S2S  N 
DMAR  =  0.0  N 

DLAR  =  0.0  N.ni  • 

Total  Impulse  -  !M0  N.s  for  t  —  0. 12  .seconds. 

EJECT2  emulates  .store  iclcn.se  from  the  inboard  station  of  the  same  Vertical  Ejector  Rack 

after  station  one  has  been  lired.  Reference  3  indicates  that  for  store  rele;i.so  from  the  E/A-  ^ 

18,  ejection  from  the  inboard  station  imparts  initial  roll  and  pitch  rates  to  the  store.  The 
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incremental  forces  and  nioinonts  are  reset  to  zero  in  the  simulation  after  0.12  seconds  have 
elapsed.  The  values  of  the  incroineiital  forces  and  moments  used  are  : 

DZAR  =  8200  N 
DMAR  =  19,8  N.m 
DLAR  =  Od.;)  N.m 

Total  Impulse  =  OSd  N.s  for  t  =  0.12  seconds. 

These  viducs  correspond  to  tiio  correct  impulse  for  release  given  in  Reference  3,  and  arc 
intended  to  provide  simple  rcprc.sentations  of  store  ejection  characteristics. 

Neither  IDJCCT  or  C.1ECT2  have  any  subordinate  subroutines. 

Subroutine  MAKEFILE 

The  MAKEFILE  snl)routiiie  is  used  to  record  the  trajectory  and  orientation  of  the  store  in 
ASCI!  rorinat  for  earli  .simulation  run.  The  sul)roiitine  uses  standard  FORTRAN  code  to  open 
(or  rni.pcn)  a  lile  at  the  lioginning  of  a  rnn,  overwriting  any  existing  contents  of  the  file  with 
the  current  trajectory  data.  Fhe  information  recorded  is  the  current  time,  x,y,r.  positions, 
yaw,  pitch,  roll  angles,  angle.s  of  attack  and  side, slip,  airspeed,  rate  of  cliange  of  altitude,  and 
normal  acceleration,  witli  the  following  format: 

'  write(l,1000(E,*  x,  y,  Paid,  Thetad,  Phid,  Mphad,  Betad, 

vde,  AN,  0,  0,  0  ", 

.v‘>crc  tlie  In.sl  three  zeroes  are  unused  tli.aniiels. 

This  Me  can  ho  used  to  obtain  the  jiosition,  speed,  and  attitude  of  the  store  in  any  post-run 
processing  that  is  re(|uired. 

The  I  ;st  time  tliat  tlie  MAKEFILE  subroutine  is  called,  the  user  is  prompted  with  the  screen 
iiiGSs.ige: 

“  At  what  time  would  you  like  to  begin  the  datafile?  "  , 

at  t  ,.ich  point  tlie  program  awaits  a  numerical  input  from  the  user.  This  is  for  the  purpose 
of  a  <  ling  ihe  store  trajectory  data  to  part  of  the  flight  path  of  an  aircraft,  so  that  the  store 
and  '  rcraft  coincide  until  the  time  specified  by  the  user  for  the  .store  simulation  to  begin. 
Thr  i;!e  is  overwritten  c.ach  lime  a  run  is  commenced,  so  that  only  the  final  simulation  run  is 
rccc  led.  The  output  file  is  called  STORE.DAT. 

Note;  For  viewing  the  recorded  trajectory,  a  graphics  replay  system  has  been 
developed  using  .th  IKLS  graphics  work  station.  The  data  file  .STORE.DAT  can  be 
transferred  to  llic  IRIS  workstation  and  tlie  trajectory  of  the  store  can  bo  viewed 
using  the  graphics  replay  program. 

Tlio  MAKEFILE  subroutine  has  no  .suiiordinatc  subroutiiie.s. 


Subroutine  SOUNDSPEED 


This  subroutine  is  (..ailed  from  the  main  ACSL  program,  with  the  current  altitude  as  input, 
and  speed  of  souikI  at  titat  altitude  as  output.  This  value  is  then  used  to  calculate  the  Mach 
number  which  is  required  by  the  AE1102  subroutine  and  the  specialised  database  extraction 
subroutines. 

•  •' 

The  SOUNDSPl'.IDU  subroutine  has  no  subordinate  subroutines. 


Subroutine  EULBRANGLES 

This  subroutine  i.s  (Nailed  from  AiCIl02  and  has  angle  of  attack  and  sid(}sli|>  as  its  inputs. 
The  outputs  iiio  the  liuler  angles  Theta  and  Phi,  with  Psi  being  set  at  zero.  The  Euler 
angles  describe  the  orioutatioii  in  space  of  a  l)ody  relative  to  an  inertial  reference  frame,  the 
Earth-Axes.  'I’he  rclcvnnt  formulao  arc  reproduced  here  (see  Refeence  4) : 


cufiO  =  cc>s(ti  )ct>.>(/7)  and 


,1111  (77) 


The  correct  sign  of  (I  is  doloi  iiiiui'd  by  testing  for  the  sign  of  n,  and  equating  the  .sign  of  9. 

Tlic  angles  of  roll  ( Phi)  and  pilcli  (Tlicta)  are  returned  to  the  AEIt02  subroutine,  wliore  they 
arc  needed  as  input  to  l>  o  database  access  routines,  along  with  the  Mach  lui tuber. 


The  .subroutine  l';iII.EllAN(ti..E.S'  has  no  subordinate  suhronliites, 


•  ♦ 


•  « 
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Subroutine  ACCES.S 

Tile  subi'oiitiiic  ACiCKS.S  is  called  from  the  beginniug  of  the  main  AC'Sl.  simulation  program. 

The  purpose  of  ACM  iES.S  is  to  provide  a  platform  for  the  calling  of  the  dalnltase  routine,  RTI),  ^ 

which  opens  and  prepaics  tlic  database  for  retrieval  of  specific  data. 

The  subroiiliiie  APCES.S  has  the  subroutine  RTI)  as  a  subordinate  .siibioutino, 
whicli  is  part  of  tlic  (laiaba.se  extraction  routines. 
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Appendix  C  -  Database  Conversion  Program 


The  database  conversion  program  was  riovolopcd  to  transform  the  format  of  the  original 
database  into  tlic  format  rc(|Uircd  by  the  standard  database  routines.  The  conversion  was 
achieved  by  reading  the  values  of  tlic  variables  from  the  original  database  which  is  stored 
as  an  ASCII  Hie  on  a  3.5  inch  ilo|)py  disk,  and  rewriting  them  to  a  new  file  in  tlie  required 
format.  Included  in  tliis  Appendix  is  a  copy  of  tlie  Database  Conversion  program. 


<8) 


program  convertformat 

real  halfl  (12) ,  half2  (5) ,  IJno  (15)  ,  13, 13,  6) 

Intogor  1,  j,  )t,  l,p,q,  r, count 


c  lump(a,a,a,  1,2,3,  4,  5,6,7,e)  refer  to  c>.cy,c'z,cl,cni, cn,olp,cmq 

print*, 'Reading  Data!' 

open  (10,  file-'m)ce2cx.dat ' ,  status-  'old',. -votlon"' road' ! 
do  20  1b1,2V04,1 
road(10,100) (half  1  (i) , i-l, 12) 
road(10,110)  (hal£2(j),  j-1,5) 

do  10  )i-l,  15, 1 

if  (X.lQ .  12)  line  (k)«»haif  1  ((t) 

10  If  ()t.gt.l2)llno(l()-half2()<-12) 

luinp(line(2) ,  lino (3) , line (4) ,  l)"llne(8) 
lump(llnet2)  ,llne(3) , lino (4) ,  2) -line (9) 
lumpdino  (2) ,  lino  (3) ,  lino  (4) ,  3)  “lino  (10) 
lump  (line  (2) ,  lino  (3) ,  lino  (4) ,  4)  -linedD 
lump ( line  (2 } ,  1 1  no  (3! ,  1  Ir.Q  ( 4 > ,  5 )  =lino  ( 12 ) 
lump (lino (2) ,  line (3) , lino (4) ,  C) -lino (13) 
lump (lino (2) ,  lino (3) ,  line (4) ,  7) -line (14) 
lump (lino (2) ,  line (3) ,  lino (4) ,  8) "line (15) 

20  continue 

open  (11,  flle<’'mX82db.dat ',  status-  'new',  action- 'write', 
scarriagocontrol" ' fortran' ) 


■At) 


• 


•  4 


•  4 


•  •  4 


print*, 'CX  roiNTS' 


write  (11, 150) 
write  ( 11, 160) 
write  (11,170) 
write  ( 1 1, 180) 
write  (11, 190) 
write  (11,200) 
write  (11,210) 
write  (11,220) 
write  (11,230) 
write  ()1,.240) 
write  (11,250) 
write  (11,160) 
wrlte(ll,270) 
wrltie(ll,280) 
write(ll,290) 

write(ll,l60) 
write(ll,330) 
write  (11, 115) 
print?, 'Writing  to  file' 

print*,  'Doing  CX. ' 
count-0 
do  30 

do  31  q-1, 13, 1 
do  32  p-1, 16, 1 

writcdl,  120)  lump <p, g,  r,  1  > 
count-count ♦ 1 
if (count . eq . 5) then 
write  ( 11 ,  115) 


•  4 


28 


•  • 


count“0 

endiC 

32  cont  iiwie 

31  contlmic 

30  continue 

wrlte<ll,320) 
print*, ’Finished  CX.' 

print*, 'Doing  CY. ' 
write  <11, 340) 
writo<U,115) 
count-0 
do  33  r-1,13,1 
do  34  q-1, 13, 1 
do  35  p«l, 16, 1 

writedl,  120)  lump(p,q,  r,  2) 
count"oount+l 
if  (count .  oq.  5)  tt»on 
writodl,  115) 
coiint-O 
otulif 

35  continuo 

34  continue 

33  continue 
write  (U,  320) 
print*, 'Finished  CY. ' 

print*,  'Doing  C?.,  ' 
writedl, 350) 
write  <U,  115) 
count"0 
do  36  i^l, 13, 1 
do  37  q“l,13,l 
do  38  p°l,16, 1 

writodl,  120)  luiflp(,  q,  r,3) 
count"count+l 
if  (count .  oq.  ft)  then 
writc(ll, 115) 
count-0 
ondif 

30  continue 

37  continue 

36  continue 
write  (11, 320) 
print*, 'Finished  CZ.' 

print*, 'Doing  Cl .  ' 
writedl, 360) 
write  (11, 115) 
count-0 
do  39  r-1, 13, 1 
do  40  q-1, 13, 1 
do  41  p-l,l6,l 

write  dl,  120)  luirp(p,q,  r,  4) 
count-count +1 
if  (count .oq. 5) then 
writodl,  115) 
count'^0 
ondit 


41 

continue 

40 

cont inue 

39 

cont inuo 

write  (11, 320) 

2'7 


print*, 'Finished  Cl. ' 

print*, 'Doing  CM.' 
write(n,370) 
wrlto<n,115) 
count-0 

do  42  r-1, 13, 1 
do  43  q-l, 13, 1 
do  44  p-1, 16, 1 

write  (11, 120) lump (p, q, r, 5) 
count-count-t  1 
if (count .oq. 5)then 
writo (11, IIS) 
couut-0 
ondif 

44  continue 

43  continue 

42  continue 

writo  (11, 320) 
print*,  'Finis)>ed  CM. ' 

print*, 'Doing  ON. ' 
wrlto(ll,300) 
writG(n,115) 
oount:«0 

v1o  45  r-l,  13,1 
do  46  q" 1 , 1 3 , 1 
do  47  p-1, 16, 1 

writo (11, 120) lump (p,q, r, 6) 
count"Count+l 
( f (count .oq.  5) then 
writo (11, 115) 
count-0 
ondif 

47  continue 

46  continue 

45  continue 
wrlto(ll,320) 
print*, Tinishod  CN. ' 

print*, 'Doing  Clp, ' 
write(ll,390) 
writo(ll,115) 
count-0 

do  46  t-1, 13, 1 
do  49  q-1, 13, 1 
do  50  p-1, 16, 1 

write (11, 120) lump (p,  q,  r, 7) 
count-count ♦! 
if (count. eq. 5) then 
wrltodl,  115) 
count-0 
ondif 

50  continue 

49  continue 

40  continue 

write (11, 320) 
print*,  'Kiivt.slied  cip. ' 

print*,  'Doing  Cmq, ' 
write  ( 1 1,  400) 
print*, 'OKI ' 
writo  (11, 115) 


# 


♦ 


•  • 


•i 


• 


•  •: 


•  • 


30 


•  • 


53 

52 

5l 


100 

110 
«  «  - 
*■10 

120 

ISO 

160 

IVO 

180 

190 

200 

210 

220 

230 

210 

250 

20o 

210 

280 

290 

300 

320 

330 

340 

3So 

360 

370 

360 

390 

40o 


print*, '0K2 ' 
count-0 
print*, 'OK3 ' 
do  51  r-1, 13, 1 
CIO  52  q-1,13, 1 
do  53  p-1,16, 1 
write  <11, 120)  luinp(p,  q,  r,  O) 
coiint’«count+l 
i£  <oount .eq.S)thQn 
wrltodl,  115) 

count-0 

ondlf 

continuo 

contimio 

ountimig 

write  <11, 320) 

print*, 'Finished  Cmq. ' 

format<415,  lfl5.4,  7£13,4) 
format  <31*,  5fl3. 4) 
forma u  i  'O' ) 

format  ( 't ' ,  lf:i3.4,  ',  ') 


format  ( ' 
format)'  **) 
tonnat  ( ' 
format  ( ' 
format  ( ' 
format  < ' 
format  ( ' 
format  ( ' 
format  ( ' 
format  < ' 
format  <  • 
format <'  *') 


tltlo  mkfl2cxciatabaso') 

real  o*,  oy,  osi,  c.l,  cm,cn,clp,  cmq' ) 
Chruput  cx  <aapHa,roij,mac)nio) ') 

thruput  cy  <aipha, roll.machno) ' ) 

thruput  cz  <alpna,toU,machno) ') 

thruput  cl  <alpha,roH,machno) '1 

tliruput  cm  <alpha, roll,mnt;Jino)  ' ) 

tlu-nput  cn  <alpho,roll,machno> ') 

thruput  dp  <alpha,roll,rnachno) ') 

thruput  cmq  <alpha, roll.machno) ' ) 


format <'  alpha  conPpt  0.0,30.0,2/') 
formate  roll  conPpt  “45.0,45,0,7.5/') 

ridl"  l‘ !  1  Sd ‘'2/  ■  r ®  ®  •  <5'  0  •  7,  0 .  C,  0 . 05,  0 . 9,  0 .  or.,  1 . 0,  1 . 05 

formate  *') 


format  < ' 
formate  cx 
foimate  Cy 
format  < '  c? 
format  < '  Cl 
formate  Cm 
formate  Cn 
formate  Clp 
format  <  ‘  Cmq 


+♦•) 

I'OINTS') 

POINTS') 

POINTS') 

POINTS') 

POINTS') 

POINTS') 

POINTS') 

POINTS') 


stop 

end 


Appendix  D  •  Example  Simulation  Run 
Hun  with  Ck'jui  Ui'loaso 

ill  tliiii  lir.sl  t'Xutiipli',  IIk<  store  is  tdonsed  itilr  the  free  stronin  without  tho  brnnfil  of  any 
ojoclioii  forr'os  or  uioiiiont<s.  Tlu'  opmilor  iupiit  is  shown  lion*  iu  lower  case,  hut  it  should  he 
noted  that  f  ile  pioj^ram  is  not  ease  sensitive,  The  eoiiiputer  echo  of  the  inputs  vite  shown  in 
upper  c.ise.  The  snurte  code  hie  is  luinied  atUu.acsl,  with  the  setup  file  natiied  adui.setup. 
'I'he  oxociitnlile  file  is  culled  iidni.exe.  ’I\i  run  the  program,  type 


1 ,  adm 

The  pioeraiii  pioinpts  wilh 

Ito.ading  dat.i . ploaso  hold 

The  prngrniii  then  .scans  the  hienkpoiiit  data,  and  rends  tlio  points  CX,  C'Y,  . 
.  .  CM(J.  When  (1M(}  has  liceii  read,  the  cursor  awaits  the  next  user  input. 


s  iliidsll)  iiisiriu'l.s  llie  piogniiii  to  take  its  iiipiil  from  the  setup  iile,  until 
ilisliurted  olherwi.se  liy  a  .sl.iteiueiit  coiilaiiied  within  tlie  input  file. 

The  piiii>  At '.SI,  seeiiioiii  of  tlio  program  pimiplirases  eiuli  instruction  on  tho 
screen,  coming  li.uk  with  S  ('M|)=l(l. 

;i,  !i  xO-lOOO 

II  llie  Miiiuhiiion  oul|iiit  file  S'l'OII IvIlAT  is  to  he  used  in  a  (light  path  r»“* 
consirnrtion  in  coiiiumliou  witli  .another  aircr.afi,  the  store  trajectory  must 
begin  ill  tlie  <lo\vM  rinige  tll.stiiiice  of  the  airciaft  at  the  time  of  release,  so  that 
the  .stole  dues  not  h'iip  liack  lo  X— (I. 

•I.  odil 

The  IfDIT  .siihroiiiino  is  not  a  siiuulard  .Af  ’SI,  coiumaiid.  hut  allows  the  sotting 
of  the  iiiici  iifi  piiriiiiielei:;  and  release  conditions.  The  aircraft  p.ar.aiiietcis  .are 
all  preset  to  values  rcliccliiij'  the  characteristics  of  the  Mk-82  store,  and  the 
release  cundilioiis  are  sol  to  default  values  which  need  to  be  modiiied  foreiirh 
ililleronl  relea.se. 

.I.  alt0=1000 

A  I,T()  is  the  initial  altitude  in  metres  of  the  .store.  N»>te  that  llie  values  within 
I'.nri'  arc  input  dillereiiily  lo  (liose  values  input  in  the  main  Ad.Sb  program. 

ti.  <10=0.01 

ot)  i^  ilir  initial  pilch  r.ile  e.sperieuced  Iry  til'' store,  in  radians  per  second. 

7.  vokO-300 


I 


■i) 


• 


•  •; 


•  # 


#  • 


•  • 


•  • 


•  •\ 


:5’.t 


•  i 

j 


•  • 


vekO  is  llic  iiiitiiil  velocity  of  llic  store  at  the  point  of  release,  in  knots. 

8.  thotd0=10 

tlietdO  is  the  initial  pitch  attitude  of  the  store,  in  degrees. 

9.  alphd0=2 

alphdO  is  the  initial  angle  of  attack  of  the  store,  in  degrees. 

1C.  betad0=2 

hetadO  is  the  initial  angle  of  .sideslip  of  the  store,  in  dcgrce.s. 

11.  trim 


The  first  time  this  routine  is  called,  the  user  is  prompted  with  the  question,  “At 
what  time  would  you  like  to  begin  the  data  file?”  The  program  awaits 
the  input  of  a  iiuitibcr  that  represents  the  time  in  seconds  that  the  piircnt  air¬ 
craft  Hies  before  the  store  is  rclea.sed.  This  is  only  relevant  if  the  flight  gr.aphic 
replay  is  going  to  he  used  in  coiijuuclion  with  another  aircraft.  If  the  store  is 
intended  to  ho  visualised  exclusively,  the  start  time  can  be  set  to  zero.  Note 
that  tliis  value  alfecls  only  the  .STORE.DAT  file  which  i.s  used  in  later  graphic 
replay,  and  not  the  simulation. 

The  only  reason  that  this  store  is  trimmed  is  so  that  the  EIGEN  function 
is  applied  to  a  liody  in  cquilibriuiii.  TRIM  introduces  artificial  incremental 
forces  and  iiiomciils  which  negate  all  accelerations  experienced  by  the  store 
at  llio  point  of  rolo.'iso, 

EDIT  has  now  been  exited,  and  the  user  is  returned  to  the  pure  ACSL  envi- 
roiiiiienl. 


•  # 


EIGEN  perforiii.s  an  eigen  analysis,  producing  eigen  vaines,  eigen  vectors,  and 
the  .l.acciliiaii  iiiatri.x.  These  may  lie  printed  out  by  inputting  the  statement ‘s 
prn=9  ’  prior  to  the  input  of  the  EIGEN  statement.  This  causes  the  output 
of  the  EIGEN  routine  to  he  diverted  to  an  output  file,  ADM.I. ,  which  may 
then  be  ])iinted  out  with  the  command  PRINT  ADM.L  after  tlie  program  has 
been  o.xitod. 

13.  dtz 

To  remove  the  incremental  forces  and  moments  quickly,  they  can  all  be  set  to 
zero  with  the  stateinciit  DTZ,  which  simply  stand.s  for  “Deltas  to  zero”. 

14.  s  tstcp*10 


The  length  of  the  siiiinlation  run  is  set  by  the  statement, 
s  tstop=10 

,  where  the  time  is  entered  in  .seconds. 


15.  go 


CO  b<>giiis  the  sitiiiilalioii. 

16.  d  alt 

Tlic  statement  ‘d  alt’  displays  the  final  value  of  the  specified  variable.  In 
thi.s  rase,  the  final  altitude  is  feet,  indicating  that  the  run  time  is  too 
short.  Repeat  the  last  three  statements  iteratively  until  satisfied  with  the 
final  altitude  of  the  store. 

When  complete,  tyi>o 

stop 

or  Coi'irol-c,  which  will  exit  the  ACSL  program.  The  STORE.DAT  file  con¬ 
tains  only  the  final  siinnlalion  run,  and  will  be  overwritten  every  time  ADM 
is  run,  .so  if  the  lilc  is  needed,  it  is  necessary  to  copy  STORE.DAT  to  a  more 
appropriate  file  nan’.e. 

If 


s  prn“0 

ha.s  been  .sol,  ADM. I.  will  contain  tlio  results  of  the  eigen  analyiis,  which  may 
then  he  printed. 


Run  Witli  Ejection  Release 

To  run  a  siinnlalion  incindiiig  simulated  ejection  forces  and  moments,  the  user  mu.st  invoke 
one  of  the  ejection  roiiiiiio.s  after  the  store  has  been  trimmed  and  all  the  incremental  forces 
set  to  zero  (step  l.'l).  Tlie  procedure  E.IECT  simulates  the  ejection  forces  and  moments 
acting  on  a  store  released  from  one  station  on  a  vertical  ejector  rack,  and  EJECT2  simulates 
ejection  forces  and  moineiits  froi-i  a  s  ^ond  station  on  the  same  rack,  after  the  first  store  is 
released.  The  applied  iinpiil.se  is  dilfcrent  in  both  eases  due  to  the  different  stiffnesses  and 
inertias  encountered  liy  the  rack  in  eacii  case,  resulting  in  different  forces  .nd  moments  being 
applied  (sec  .Appendix  It. '2  for  details  of  ejection  impulses).  The  actual  ojection  of  the  store  is 
completed  in  0.12  seconds,  so  the  ACSL  program  tests  for  time  greater  than  0.12  seconds,  and 
sets  '.he  increiiiciilal  forces  to  zero  wlicii  this  condition  is  true.  Dccansc  of  this,  it  is  necessary 
to  call  EJECT  or  E.1ECT2  each  time  the  simulation  is  run,  ie  before  the  command  GO  is 
input. 

Note;  if  it  is  necessary  to  trim  the  store  again  during  the  course  of  a  simulation  study,  the 
program  will  crash  and  c.xit  unless  at  least  one  of  the  incremental  forces  or  moments  is  rc.set 
to  a  non  zero  value  to  prevent  a  division  by  zero. 
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Mk>82  bombs  CATcoofUGs 

Tret  flight  1902 

Flight  dynamics 
Exlemal  store  senaration 


H.  AnriMcr 

.4  six  degree  of  freedom  flight  dynamic  model  incorporating  an  aerodynamic  database  derived  from  DSJO 
wind  tunnel  data  was  developed  for  the  Mk-82  store.  The  model  can  be  linked  with  aircraft  flight  dynamic 
models,  and  forms  part  of  a  general  programme  to  model  aircraft  carriage  and  delivery,  weapon  release 
and  store  ballistics.  A  rudimentary  ejection  model  was  also  incorporated  to  represent  the  forces  and 
moments  experienced  by  the  store  at  the  point  of  release.  The  model  is  valid  for  free  stream  flight  only,  and 
ignores  the  effects  of  aircraft- store  inieroclion. 
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